Three-fermion problems in optical lattices 
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We present exact results for the spectra of three fermionic atoms in a single well of an optical 
lattice. For the three lowest hyperfine states of 6 Li atoms, we find a Borromean state across the 
region of the distinct pairwise Feshbach resonances. For 40 K atoms, nearby Feshbach resonances are 
known for two of the pairs, and a bound three-body state develops towards the positive scattering- 
length side. In addition, we study the sensitivity of our results to atomic details. The predicted 
few-body phenomena can be realized in optical lattices in the limit of low tunneling. 
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Introduction.- Experiments with cold atomic gases 
make it possible to study strong-interaction physics in a 
controlled manner. When an atomic gas is loaded into an 
optical lattice, typically a few atoms reside in each well. 
Therefore, optical lattices can be used to investigate few- 
body phenomena when the tunneling barrier between po- 
tential wells is high 1]. In dilute gases, the interactions 
are governed by the S-wave scattering length a, which 
can be tuned across atomic Feshbach resonances. Con- 
sequently, three-fermion problems in optical lattices can 
access nearly the entire landscape of fascinating few-body 
phenomena. When all scattering lengths are large, the 
few-body physics of dilute gases exhibits universal prop- 
erties because there are no length scales associated with 
the interaction. These universal aspects stretch across 
physics: For example, the large scattering-length physics 
predicts a linear correlation for ground-state energies of 
few- helium clusters or light nuclei Q. 

The first step towards realizing isolated few-atom sys- 
tems was the formation of molecules from fermionic 
atoms in an optical lattice In addition, the ground- 
state energy of two particles in a single well of an op- 
tical lattice was measured across a Feshbach resonance 
by rf dissociation to noninteracting fragments. The 
measured ground-state energies (for varying scattering 
length) agree very well with the theoretical prediction 
for two particles in a harmonic oscillator potential Q . 

For three identical bosons, the scattering lengths of all 
pairs can be tuned with one Feshbach resonance and the 
three-boson problem is universal: On resonance, there is 
an infinite tower of Efimov trimer states with consecutive 
binding energies E n /E n+ i « 515 in free space Efi- 
mov states become bound for a finite, negative scattering 
length, and lead to resonances in the rate for three-body 
recombination. Recently, the first Efimov resonance was 
observed in a gas of cold Caesium bosons [5(. For three 
particles in a harmonic well, the spatial confinement re- 
stricts the size of the Efimov states, and thus the accu- 
mulation of bound states does not persist. However, Stoll 
and Kohler have shown that it is possible to study the 
first Efimov state isolated in a single well of an optical 



lattice [oj . This three-boson problem was also studied by 
Jonsell et al. in an adiabatic approximation 0]. Efimov 
trimcrs can be examples of Borromean systems (when the 
scattering length is negative): three-body bound states 
for which all pairs are unbound. Other Borromean sys- 
tems are the 6 He nucleus and the 11 Li halo nucleus ( 4 He 
or 9 Li and two neutrons) Q. 

For fermions, in contrast to identical bosons, the Pauli 
principle restricts S-wave interactions to different hyper- 
fine states, which have distinct pairwise Feshbach reso- 
nances. In this Letter, we calculate exactly the spectra 
of the three lowest hyperfine states of 6 Li or 40 K atoms 
confined to a single harmonic well of an optical lattice, 
using effective field-theory interactions. For 6 Li atoms, 
there are Feshbach resonances between all pairs 9] , two of 
which overlap closely, and we find a Borromean state that 
extends across all Feshbach resonances. For 40 K atoms, 
Feshbach resonances exist between two pairs, whereas the 
third pair interacts weakly. In this case, we find that 
a three-body state only develops towards the positive 
scattering-length side, and is bound around the doubly- 
interacting particle. The spectra for three 6 Li or 40 K 
atoms in optical lattices are very rich, and the Bloch- 
Horowitz method employed here is ideally suited to iden- 
tify the angular momenta of the states. The predicted 
few-body phenomena can be realized in optical lattices 
in the limit of low tunneling. 

Three-fermion problems. - Across a Feshbach resonance 
the dependence of the scattering length on the mag- 
netic field B is given by a{B) = <z bg (l - A/[B - ~B\) 
where B and A are the position and width of the reso- 
nance and ab g denotes the background scattering length. 
For 6 Li, the three trapped hyperfine states are the low- 
est magnetic sub-states: |1) = \F, rap) = |l/2, 1/2), 
|2) = |l/2, —1/2) and |3) = |3/2,-3/2), with distinct 
Feshbach resonances, as shown in Fig. Q] The resonance 
positions are B 12 = 83.41 mT, B 13 = 69.04 mT and 
B23 = 81.12 mT Q, and we use the Feshbach parame- 
ters plus leading-order correction determined in [9j . The 
relevant hyperfine states for 40 K are |1) = |F, mj?) = 
19/2,-9/2), |2) = |9/2,-7/2) and |3) = |9/2,-5/2). 
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As shown in Fig. [2j nearby Feshbach resonances are 
present between the states 12 and 13, with Feshbach 
parameters: B n = 202.10 G Jj$, A 12 = 7.8 G [3, 
B 13 = 224.21 G, A 13 = 9.7 G [llfand background scat- 
tering length dbg ~ 174 ao [HI ]. 

We work with effective field-theory contact interactions 
regulated by separable cutoff functions [l3| 

V( P ' 1 p) = —g(B 1 A)e-« 2 +^\ (1) 
m 

where p and p' denote incoming and outgoing rela- 
tive momenta, and the coupling g(B,A) is determined 
from the scattering length through g(B,A) = a(B)/(l — 
Aa(B) I V2tt) . If a is weak (a ~ R, where R is the range 
of the interaction), it is possible to choose the cutoff in a 
wide range with |A<x| -C 1, and one recovers the standard 
pseudo-potential for low momenta V(0,0) = 4Trh 2 a/m. 
The cutoff generates an effective range r c ~ 1/A and 
higher-order terms, which we render small with large cut- 
offs. In addition we can vary A. This probes neglected 
effective range effects and sensitivity to atomic details. 

For the separable interaction, Eq. (fl|), the two-body 
problem in a harmonic oscillator potential can be solved 
exactly. Following Busch et al. [14J, the intrinsic energy 
E = e fku is given by 

171 /3 3 c.7 e . [ l-x 2 1 2 \ [7T- 
2-^2' 4 ~ 2> 4 ~ 2> Ll+^J ) , V27T 
5 V x = TT i 2) 

(i + ^) 3 (!-!) a / fo 

where 2-F1 is the hypergeometric function, x — Ab and 
b = ^Jh/moj is the oscillator length. For large cut- 
offs x — > 00, we recover the result of Busch et al, 
\/2T(§ - |)/r(i - § ) = 6/0 Q. The energy from 
Eq. © is within 3% (or 7%) of the latter for Ab = 100 
and all scattering lengths except the tight-binding region 
< a/b ^ 1 (or 0.5). Typical well frequencies in opti- 
cal lattices are v ~ 100 kHz. Consequently, the oscillator 
length b ~ 1000 ao (a denotes the Bohr radius) is large 
compared to the range of atomic interactions R ~ 10 ao. 
Finally, the measured ground-state energies of two par- 
ticles in a single well of an optical lattice agree very well 
with this result [3j even down to \a/b\ ~ 1. 

Block- Horowitz method. - The spectrum for the intrin- 
sic energy of three hyperfine states in a harmonic well is 
determined from the three-body Hamiltonian 

3 

H = H + V = J2 H (ri) - ff ,cm(R) 
i=l 

+ Vi 2 (ri,r 2 ) + Vi 3 (ri,r 3 ) + V^3(r 2 ,r 3 ) , (3) 

where the noninteracting part is given by Hq(t) = 
-fi 2 V 2 /2m + muj 2 r 2 /2 and R = (r x + r 2 + r 3 )/3 de- 
notes the center-of-mass (cm) coordinate (J?o,cm uses 
mass 3m). Within Jacobi coordinates, the eigenstates 
of Hq are characterized by 

#o|12,3) = (JV + 3)/ja>|12,3). (4) 




B [mT] 

FIG. 1: Top: Feshbach resonances between the three lowest 
hyperfine states of 6 Li atoms Bottom: The spectrum 
for the above 6 Li states in a single well of an optical lattice 
with v = 270 kHz versus magnetic field. These BH results are 
for A& = 100. The three-body energy E3 is measured from 
the three-body dissociation threshold given by the minimal 
two-body energy i?2,mm (from Eq. (J3J) of the pairs and the 
additional particle in the noninteracting ground state E\ = 
3/2/kj (5/2/wj for the negative parity state P=-l). For B > 
(<)73.2mT, i52,min is determined from 012(023). The state 
labeled B is the Borromean state, and En_L denotes the n-th 
excited state with angular momentum L. 



Here, 12 = ni2, ^12, mi2 are the radial, angular and mag- 
netic quantum numbers of the pair 12, and 3 = 713, I3, 
refer to the quantum numbers of the third particle with 
respect to the cm of pair 12. We classify these noninter- 
acting states \Ni) according to their principal quantum 
number N = 2n\2 + l\2 + 2ri3 + Z3, where the subindex i 
denotes all possible states for fixed total N. 

We solve the three-body problem (Ho + V)\ip(E)} = 
EW)(E)) using the Bloch-Horowitz (BH) approach [l5|, 
Il6j . The BH method diagonalizes an effective Hamil- 
tonian in a truncated space of P = J2N i <N ln 

\Ni)(Ni\ 

low-energy excitations, such that the low-lying spectrum 
is exactly reproduced. Inserting 1 = P + Q, we obtain 
the projections of the three-body Schrodinger equation: 

P(H + V)(P + Q)\iP(E))=EP\iP(E)), (5) 
Q(H + V)(P + Q)\ij(E)) = E Q\HE)) . (6) 

Since [P, H ] = [Q, H ] = and PQ = 0, we can solve 
Eq. © for Q\ip(E)) ={E- H Q )- l QV\^{E)) and insert 
the latter into Eq. (O . This leads to the equivalent prob- 
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FIG. 2: Top: Feshbach resonances between the three lowest 
hyperfine states of 40 K atoms [10-12]. Bottom: The spectrum 
for the above 40 K states in a single well of an optical lattice 
with v — 100 kHz versus magnetic field. These BH results 
are for a23 = and Ab — 100. The three-body energy Es is 
measured from the three-body dissociation threshold given by 
the ground-state energy of pair 13 (£2,(13)) an d the additional 
particle in E\ = 3/2hu> (5/2huj for P=-l). The state labeled 
G is the three-body ground state and EnL is as in Fig. [TJ 



lem in the truncated space, 



P[Ho + V + V 
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V)P\^{E))=EP\i>{E)), (7) 



where the effective Hamiltonian H c g(E) (given by the op- 
erator in parentheses in Eq. (7])) depends self-consistently 
on the exact energy and exactly reproduces the low-lying 
spectrum, as long as the eigenstate has overlap with the 
truncated space P\ip(E)) ^ 0. Finally, we use a Faddeev 
decomposition \i/>(E)) = (1 + P12P13 + Pi 2 P23)|V>(P))i2 
to construct H e g(E) (P^- is the permutation operator). 

The good quantum numbers of the interacting eigen- 
states are parity P, total angular momentum L and pro- 
jection L z : \ip{E)) = \E; P, L, L z ). Here we solve the BH 
Eq. |7]) in an uncoupled basis for L z = (mu = m 3 = 0). 
Since \E;P,L Z ) = £ L C L (P; P, L z ) \E; P, L, L z ) , the re- 
sulting spectra automatically contain states with all pos- 
sible angular momentum quantum numbers. We will 
use the BH overlap condition to identify their angular 
momenta. The BH method has been used to calculate 
the ground-state properties of light nuclei [3] , and as a 
check, we have reproduced the results of Stoll and Kohler 
for three identical bosons Q. For the separable interac- 
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FIG. 3: The cutoff dependence of the Borromean and the first 
excited state of three 6 Li atoms versus magnetic field. 



tion, Eq. |T]), it is possible to calculate the necessary BH 
two-body matrix elements analytically. In addition, we 
have found it sufficient to keep Z 12 , Z3 ^ 3. 

Results.- The magnetic field dependence of the spec- 
trum for three 6 Li atoms in an optical lattice with v = 
270 kHz is shown in Fig. [TJ The BH results are indepen- 
dent of iV max . In particular, all states are present in the 
lowest 7V max = calculation (with Z12 = £3 = 0), and thus 
all positive parity states shown have angular momentum 
L = 0. The lowest negative parity state has L = 1. 

We find a deeply-bound Borromean state B that exists 
on the negative scattering length side and extends across 
the Feshbach resonances. Note that there are many very 
deeply-bound two-body states present. This state can be 
viewed as a collective state within a schematic model [13] ■ 
For B > 75 mT, the excited and negative parity states 
depend very weakly on the magnetic field in Fig. [TJ since 
the Feshbach resonances of a\i and 023 are very close and 
here this two-body energy is subtracted. The first excited 
state Elo is adiabatically connected to the noninteracting 
N = state at high magnetic fields. Similarly, the states 
E2o and E3o connect to the two N — 2 states with Z12 = 
£3 = 0, where the other three states of the noninteracting 
N = 2, L z = multiplet are higher in energy since they 
are less sensitive to S-wave interactions. 

In Fig. [2] we show the spectrum for three 40 K atoms 
in an optical lattice with v = 100 kHz versus magnetic 
field. In this case, there are two Feshbach resonances be- 
tween pairs 12 and 13, and we have taken the third pair 
to be noninteracting 023 = in this calculation. We find 
that a bound three-body state only develops towards the 
positive scattering-length side of both resonances. The 
three-body state is bound by the doubly-interacting par- 
ticle 1. The qualitative features of the spectrum do not 
depend on 023 for | CE.23 / ao | ^5 100. For instance, with a 



4 



a 23 -174 a , Ab=100_ _. 




Ab=600 

Ab=400 

b=950a n Ab=200 

v=100kHz Ab=100 



210 215 
B[G] 



220 225 230 



FIG. 4: The cutoff dependence of the ground state of three 
40 K atoms versus magnetic field. 



repulsive 023 = ab g ~ 174 ao, the spectrum is moved up 
and the ground state becomes bound only for lower mag- 
netic fields (see Fig. 0] note that the precise value of 023 
has not been calculated). 

Two of the states of Fig. [2] (E2 2 and E4 2 ) are not 
present in a iV max = calculation, but exist for all larger 
JV max ^ 2, and thus have angular momentum L = 2. At 
high magnetic fields, we recover the five states adiabati- 
cally connected to N — 2 (note Z13 = 1^ = 1 can couple to 
L = 0, and also the avoided level crossing for two of the 
L = states). For B = 209.9 G, we have ai 2 (B ) = 0, 
and the only interaction is for pair 13. Therefore, the 
low-lying states are E3 — ^2,(13) — 3/2huj = 2n 2 + l 2 for 
the positive parity state (the excitation of the cluster 13 
comes higher in energy; and 2ri2 + I2 — 1 with —b/2fno for 
P=-l), in agreement with Fig.fJ The states E2 2 (l 2 = 2) 
and P=-l (l 2 — 1) follow only the two-body energy of 
the right Feshbach resonance. Finally for B > Bo, the 
interaction between 12 becomes weakly repulsive, which 

„, 40 



requires Ai 2 a < \/2tt. Our °K results are for Kb = 100, 



except Ai20b g = 1 for B > Bo. For Ai2<2 < v27r, we find 
a very weak cutoff dependence from the repulsive part of 
the 12 interaction. 

We can vary the cutoff and thus probe the depen- 
dence of our results to the effects of an effective range 
and many-body interactions [l3| . In Figs. [3] and |3J we 
show the cutoff dependence of the Borromean and first 
excited states for 6 Li, and the ground state of the 40 K 
three-body problem. While this excited state (and all 
others) are well converged, the ground state energies con- 
verge slower and show a sizeable dependence on A6, as the 
binding energy increases. Therefore, these states become 
sensitive to further atomic details, such as the effective 
range. Our results also indicate that there is no limit 
cycle in a harmonic oscillator potential with b/R < 100, 
and thus the power-counting of three-body interactions 
in the corresponding effective field theory (EFT) must 
change compared to free space. This may be important 
for a pionless EFT [l3[ for nuclei in an oscillator basis. 

Conclusions.- Optical lattices open a frontier to con- 



trolled strong- interaction few-body physics, in addition 
to simulating condensed matter models. In this Letter, 
we investigated three-fermion problems in optical lat- 
tices for 6 Li and 40 K atoms using the BH method. For 
6 Li atoms, we find a Borromean state on the negative 
scattering-length side that extends across the Feshbach 
resonances. In contrast, for 40 K atoms, one of the pairs 
interacts non-resonantly at the relevant magnetic fields, 
and a three-body state, bound by the doubly-interacting 
particle, develops towards the positive scattering-length 
side. While the quantitative results of the ground states 
are somewhat sensitive to atomic details, these features 
are independent thereof and also independent of the pre- 
cise oscillator frequency. The three-fermion spectra are 
very rich, and we have identified the nature and angu- 
lar momenta of all low-lying states. We predict a Bor- 
romean state in optical lattices under the conditions of 
overlapping or close Feshbach resonances for all pairs and 
attractive background scattering lengths. 
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